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1. Introduction 

It has been argued that the QCD dense matter in relativistic heavy ion collisions (HIC) 
exhibits properties of a collective fluid-like motion with low viscosity to entropy ratio (see 
e.g. Q). Thus, relativistic hydrodynamics has become an important analysis tool for HIC. 
Relativistic hydrodynamics is formulated in terms of conservation laws of the stress-energy 
tensor and various conserved currents. 

It has been recently revealed that the hydrodynamics description exhibits an interesting 
effect when a global symmetry current of the microscopic theory is anomalous. This has 
been first discovered in the context of the gauge/gravity correspondence j|, Q. The 
Chern-Simons term in the gravity action, which corresponds to having an anomalous global 
symmetry current in the dual gauge theory, has been shown to modify the hydrodynamic 
current by a term proportional to the vorticity of the fluid. 

At first sight the additional vorticity term seemed to contradict the second law of 
thermodynamics ||. This, however, has been resolved by a redefinition of the entropy 
current in ||. In this work we suggest an experimental signal, which is a consequence of 
the anomaly effect. We will consider the effect of the vorticity term as well as that of the 
gauge fields. 

The major effort in the experimental study of QCD topological effects in the context 
of HIC has been focused on charge separation. The origin of this effect is the assumption 
that in the deconfined phase of QCD a non-trivial, space-dependent, value for the QCD 
angle can be generated. In this P violating vacuum a strong magnetic field would induce an 



- 1 - 



electromagnetic current along the magnetic field lines. The experimental signature of this 
effect is an asymmetry in the charge distribution of the scattered particles in non-central 
collisions (?|, ||. 

In this paper we propose an observable which is charge independent. The basic idea is 
that the axial charge density, in a locally uniform flow of massless fermions, is a measure of 
the alignment between the fermion spins. When the QCD fluid freezes out and the quarks 
bind to form hadrons, aligned spins result in spin-excited hadrons. The ratio between 
spin-excited and low spin hadron production and its angular distribution may therefore 
be used as a measurement of the axial charge distribution. Due to the short lifetime of 
high-spin hadrons such as the p mesons and A baryons, we propose to focus on narrow 
resonances such as . We will predict the qualitative angular distribution and centrality 
dependence of the axial charge. Our main proposal is that for off-central collisions we 
expect enhancement of Q~ production along the rotation axis of the collision (see fig. [l]). 




Figure 1: The left figure shows an off-central collision of two Gold ions (the beam direction is 
transverse to the plain of the plot). In blue we see the two spatial angles in which we compare the 
production rates - in the upper and lower " cones" we expect to find an enhancement of spin excited 
hadrons due to the non-zero axial charge Qa in the fluid, and the "belt" can be used to measure the 
production without axial charge. The plot on the right shows qualitatively the predicted centrality 
dependence of the effect. 

As we will explain, the calculation of the precise magnitude of the effect requires a 
detailed numerical analysis as well as making certain assumptions about the hadronization 
process. In particular, the Bjorken flow ansatz which is very useful in the numerical 
analysis of the hydrodynamics equations cannot be used in this case. Instead, we will use 
an estimate for the axial charge distribution at early stages in the evolution of the system. 

The paper is organized as follows: In section || we briefly review the theoretical back- 
ground for the anomaly effect in the hydrodynamics framework and discuss the issues 
involved in detecting this effect in heavy ion collisions. In section ||| we present an experi- 
mental signal, use the Glauber model as the initial condition to estimate the axial density 
and give an outline for the data analysis. The last section is devoted to a discussion and 
outlook. 
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2. Hydrodynamics, Triangle Anomalies and HIC 



2.1 Relativistic Hydrodynamics With Anomalous Currents 

The hydro-dynamic description of a classical relativistic fluid is a set of conservation equa- 
tions of the stress-energy tensor and the global symmetry currents 

a^ = o, d,f a = o , (2.i) 

where 

TV = (e + P)u^u u - g^P + 
f a =p a u^ + ^. (2.2) 

is the fluid velocity field, normalized such that u^u^ = —1, e, P and p a are the energy 
density, pressure and charge densities, respectively, r^, Va are the dissipative terms that 
contain derivatives of the various fields. There exists an ambiguity in the definition of the 
fields, which we will fix by choosing the Landau frame in which the velocity represents the 
rest frame of the energy density 1 . The conditions for the Landau frame are 

u^ v = = • (2-3) 

In the Landau frame the global symmetry current takes the form 

f a = p a u» - a a T{sT + ^u v )d v + a a EH , (2.4) 



where T, p a and a a are the temperature, chemical potentials and the conductivities of the 
medium, and E% = F^ 1 ' u v is an external field which is coupled to the current ja- This form 
of the current is modified when the global symmetry current corresponds to an anomalous 
current in the microscopic theory [Q, [|, ||, ^]. In the latter case the current takes the form 

f a = p a u» - a a T(jT + u»u»)d v (^) + a a E% + ^ + , (2.5) 



where the vorticity and the magnetic field are defined by 



B£ = \e^u v {F b ) Xp . (2.6) 
The vorticity and magnetic field coefficients for abelian currents read || 

t _ n 2 n PbPcPd 

C,a — ^abc^bPc ~ T,Pa^bcd ; — 

3 e + P 

tab = C abc p c - - PaCbcd • ( 2 -7) 



1 Another useful choice is the Eckart frame, where the velocity is determined by the rest frame of the 
charge. In this case v£ = 0. 
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C a b c is the coefficient of the triangle anomaly of the currents jajfj and j£ , 

EiQIQiQl 



C a bc 



« = -~C abc e^F b F c 



2vr 2 

flV^ Op 



and Q a is the charge of the i'th Dirac fermion with respect to the a symmetry. Note that 
we absorbed the coupling constant in the definition of the gauge fields. 

2.2 Heavy Ion Collisions and the Axial Anomaly 

A set-up in which the hydrodynamics description seems useful is the description of high 
density QCD matter created in heavy ion collisions. In very energetic collisions the hot 
dense QCD matter can go through a phase transition into a deconfined phase described 
by a fluid-like collective motion of quarks and gluons. The comparison of the relativistic 
hydrodynamics simulations with the data suggests that the relativistic fluid is characterized 
by a low shear viscosity to entropy ratio, which is a property of strongly coupled systems. 
In this work we will relate to experimental observables the effect of triangle anomalies on 
the hydrodynamics description discussed in the previous section. 

We will consider a deconfined QCD fluid phase, with three light flavors and chiral 
symmetry restoration. The global U(l) currents correspond to U(1)b, the Cartan subal- 
gebra of SU(3) of flavor (which we denote by f/(l)/3 and U(l)s), and their axial versions. 
Accordingly, the relevant currents will be denoted by and , jj 5 , j^ 5 . The 

electromagnetic current will be considered as a linear combination of the vector currents, 
J7 = ij + \{3b j's)' ^ n ^ n ^ s wor k we are interested in the axial current for which the 
relevant triangle anomalies are 

13 3 11 

Cabb = Caii = ^j, Cass = Cabs = C All = ^ ( 2 - 9 ) 

The vorticity coefficient will therefore be given by 

U = ^2^1 + \ti + 34 - + 0{p A ) (2.10) 



where, as will be discussed in section |3.3| , we neglect the terms in the vorticity coefficients 
proportional to the axial density because they are subleading. The only external magnetic 
field coupled to these currents is the electromagnetic field, therefore we will use the linear 
combination of these global currents which couple to this field 

d 7 = ^7 + 0(p A ) = ^2 (2/x/ + PB + US) + 0{ PA ) (2.11) 

In the deconfined phase there are in principle additional degree of freedom that must 
be taken into account, and these are the color current and gluon fields. We consider these 
as external sources and their contribution to the equations of motion for the axial current 
is in the (non-)conservation equation for the axial current 

9,3a = -\c All e^F, v F ap - \c Acc e^PGl u G a ap . (2.12) 
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where G^ v is the gluon field, and 

3 

°ACC = ^2 • ( 2 - 13 ) 

We will estimate this effect, but will not assign a nonzero chemical potential to the color 
current. 

Note that equation ( 2.12j ) requires some clarification. In a hydrodynamics description 



one considers conserved currents, while the axial current is not conserved. One can allow 
a non-conservation of the currents due to external sources in a hydrodynamics framework. 
However, can we actually consider the electromagnetic and gluon gauge fields as external 
sources? The fact that fields generated by the fluid are considered here as external, is 
similar to the standard discussion of stellar magneto-hydrodynamics, where the magnetic 
fields generated by the rotation of the star are considered as an external force. Moreover, 
in an off-central collision the gauge fields are generated mostly by spectator nucleons, and 
can be treated as external force terms. 

When attempting to make precise quantitative predictions one encounters several is- 
sues. First, we need to estimate the chemical potentials appearing in the formula for the 
coefficients £4 ( |2.1C ) and (2.11). We will provide such an estimate in 3.3, In general, in the 



analysis of HIC in the hydrodynamic framework one assumes a small chemical potential 
compared to the temperature, which can imply an effect too small to be detected. Second, 
in most analyses of HIC, one assumes a Bjorken flow ansatz, where all the observables are 
boost invariant. This means that the divergences of the vorticity and magnetic fields are 
zero, and the equations for the axial density are trivial (see eq. ( [2.1| ) and ( |2.5D ). Another 
way to see this is to note that the Bjorken ansatz assumes zero angular momentum in 
the reaction plane, while the existence of large angular momentum is the main source for 
vorticity and magnetic field in the fluid. Solving the hydrodynamics equations without 
the Bjorken ansatz is a difficult task, which we will leave for future work. Instead we will 
estimate the axial density distribution at an early stage of the flow without fully solving the 
complete set of equations. Third, there is no direct experimental access to the axial charge 
distribution of the fluid, since all the information regarding chirality is erased during the 
hadronization process. We will propose a solution to this problem in the next section. 



3. An Experimental Signal 

3.1 The Axial Charge and Enhanced Production of High-Spin Hadrons 

The phase transition from a fluid state of QCD matter into hadron gas is arguably the 
least understood stage in the hydrodynamic description of the collision. It is unclear where 
and when the phase transition occurs and how exactly the free quarks bind and form 
hadrons. Nevertheless, using a phenomenological description of the process we argue that 
non-zero axial charge in the context of heavy ion collision can lead to an enhancement in 
the production of spin-excited hadrons. 

For simplicity, we will assume that the momenta of the fermions in a small volume 
of moving fluid are pointing in the same direction (see figure ||). Note also, that in the 
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zero mass limit, a non-zero axial charge means a preferred helicity for these fermions. The 
combination of these two statements means that a non-zero axial charge enhances the 
probability for quark spins to be aligned. When the fluid freezes-out and particles with 
aligned spins bind to form hadrons, the bound states cannot have low intrinsic spin 2 . We 
therefore propose that a non-zero axial charge enhances the production rate of spin-excited 
hadrons. 

In order to make this more quantitative, let us define 

where iVj(O) is the number of hadrons of species i detected in a solid angle and Nt t(ty 
is the total number of detected particles in that angle. In order to find the dependence of 
this quantity on the axial charge in a volume element, Qa = J dVpA = Ni — Nr, we define 

} - jffgj) A *0 - jffgj) 

where N* ot (£l) is the total number of spin-excited hadrons detected in f2, and the index 
means that this quantity is evaluated at zero axial charge in the QCD fluid. Note that 
assuming a radial flow, the relevant volume of fluid is a cone covered by the angle O (See 
fig. [2]). A^ and A* are parameters which depend on the hadronization process. They are 
unknown theoretically and will be measured in scattering angles which cover cones with 
no axial charge, 

A° = A 4 (O (Qa=0) ) . (3.3) 

We now divide the total number of fermions denoted by ntot (both left handed and 
right handed quarks and anti-quark), in a given volume of QCD fluid into two groups: 
one with an equal number of left and right fermions, and the other with left handed only 
(or right handed only, if the axial charge is negative). Thus, their relative portions in a 
given volume are Utot ~\ ( ^ A \ anc [ \Qa1 respectively. While the species and spin of hadrons 
produced from the first group will be distributed according to the "regular" ratios A^ as 
dictated by the hadronization process, the second group can only bind into spin-excited 
hadrons, because the spins are aligned. The species of the particles in the second group, 
therefore, will be determined by A* . If we take, for example, the proton as a representative 
of the low spin hadrons and the A resonance as a representative of the spin-excited hadrons 
we can write 

_ ntot(Sl) - \Q A m , 

Aa(S2) = 7o\ A A H 77vT A A ( 3 - 4 ) 

where the notations ritatiP) and Qa(Q) mean the fermion number and axial charge in the 
volume covered by the angle £1. 



2 for example, a pseudo-scalar meson can only be formed by quarks with anti-aligned spins 
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spin-excited hadrons 



Figure 2: The figure on the left shows the spatial angle ft and the cone of QCD matter which 
flows in its direction. On the right we show a zoomed-in cartoon of a small volume of fluid with 
non-zero axial charge. The small circles represent the fermions and the thick arrows represent their 
spins. The larger circles represent the bound states. Bound states with aligned spins can form only 
spin excited hadrons. 



Using this result we get the following ratio: 



Ap(fi) Xp \ ntot{n)-\Q A (n)\ X° A 

A* i 

Treating as some unknown parameter of order 1 , and writing in terms of experimental 
quantities, we get 

N A (n)/N A (n iQA=0) ) \Q A (n)\ _ \Q A (n)\ \q a (q)\ 

N p (n)/N p (n (QA=0) ) °" n tot ^l) - \Q A (P)\ ~ n tot (Q) * N part (Q) [ ' ' 

In the last step we assumed that ntot, the number of particles in the fluid, is proportional 
to the number of nucleons participating in the collision, N par t, which will be discussed in 
more detail in the next section. A more precise analysis of the effect, which requires a 
numerical solution of the hydrodynamics equations, is suggested in the discussion. 

3.2 The Glauber Model and the Distribution of Axial Charge at Early Stages 

In this section we describe a method of estimating the axial charge distribution at early 
stages of the evolution of the system. We focus on the QED contribution, and relate it to 
QCD effects in the next section. The non-conservation equation for the axial current in 
the Landau frame is 

d P f A = dnipAu" - a A T(g^ + u»u v )d v (^) + + tl^) = C All E^ , (3.7) 



•-I A* JV° m 
-jft- = jfifr can naively be approximated as -y , using the degeneracy of the various spin multiplets. 



'tot 
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where the coefficients £a and are taken from (|2. 10 ) and fl2.11| ), and since there is no 
external field associated with the axial charge, we dropped the E*i term. E^ and B^ are 
the electromagnetic fields generated in the HIC. 

Under the assumption of uniformly distributed chemical potentials and anomaly coef- 
ficients d^HA^o ~ d^A\t Q ~ 0, the equation takes the form 



(3.8) 



Thus, the RHS can be considered as source terms for the "classical" axial current, pau^ ' ■ 
At t = to, the axial density is zero and so is its spatial derivative (pa = diPA = 0) and we 
get 



dtPA 



1 



to 



w 



to 



(3.9) 



Given the time derivative of the axial density we can estimate its distribution at early 
stages 



PA\t +At ~ AtdtPa 



to 



At 

u 



t (c a ^e^ - - d^Bnk , 



(3.10) 



where At is a time interval in which this linear approximation is assumed to be valid. 

The important assumptions that have been made so far are the uniform distribution 
of the chemical potentials and anomaly coefficients, and that for sufficiently short times, 
the evolution in time of the system can be approximated by a linear expansion. We will 
now estimate the source terms at the initial conditions using the Glauber model. In this 
model the energy and velocities at the initial conditions will be obtained assuming that 
the nucleus density is given by the Woods-Saxon distribution 



p(r) 



Po 



1 -)_ e -(r-Ro)/ao 



(3.11) 



where p{r) is the nucleon density, and for gold ions we use the values oo = 0.54 fm and 
Ro = 6.4 fin Q. po is determined by the condition J dVp = A = 197. It is useful to define 
the "Thickness function" 



/oo 
dzp(\/ x 2 + y 2 + z 2 ) . 
-oo 



(3.12) 



In a non-central collisions, we set the origin of our coordinate system between the centers 
of the two colliding ions, and set the y axis along the rotation axis. The number of 
participating nucleons, N part , is given by 



N P art{b) = / dxdy 



T(x+-) [ 1 



+ T(x- -) I 1 



a NN T(x 




(3.13) 
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where the nucleon-nucleon scattering cross section ctattv ~ 40 mb. The energy density of 
the fluid in the Glauber model is proportional to the product of the thickness functions 

e(x,y)\ t0 cxT(x-b/2)T(x + b/2) . (3.14) 

The initial velocity of an infinitesimal area of fluid v z (x, y) is assumed to be the center 
of mass velocity at that location 

v z {x,y)\ t0 ~ T (x + b/2,y)+T(x-b/2,y) ^ 15j 

/3 ~ 1 is the velocity of the colliding ions. The z-component of the 4- velocity, u z , is given 
by jv z . 

We now use this model to evaluate the source terms in (|3.8f). We have 



= \e^{d^u v ){d p u a ) . (3.16) 

Since at t = to all the velocities are along the z axis, namely in the beam direction, we find 
that u x \t = u y \t = 0, and ut\t = y/T+~v%. Therefore, at this stage the only non-zero 
contribution to the divergence of the vorticity is 

d^\ to = [{d t u x )(d y u z ) - {d t Uy)(d x u z )] | to . (3.17) 

In order to estimate dtu x we use the hydrodynamics equations for a perfect fluid: 

u u d u u^ = ^ ^ (3.18) 



using the equation of state for conformal hydrodynamics e = 3P (this condition can be 
easily relaxed), and the initial conditions u x = u y = d z u x = 0, we are left with 

d t u x \ to = - 1 2 -rl*o • ( 3 - 19 ) 

Similarly, we can find the divergence of the 4-dimensional magnetic field 

d^\ t0 = l -e^(d,u v )F ap \ t0 

F xz {d y u t - d t u y ) - F yz (d x u t - d t u x ) - F xt d y u z + F yt d x u z ^j \ tQ . (3.20) 

The electromagnetic fields generated by the colliding ions can be obtained by boosting the 
electromagnetic field from the ion rest frame. The electric field can be estimated assuming 
the charge density is also given by the Woods-Saxon distribution. 

3.3 Estimating the Magnitude of the Effect 

There are several additional factors that we have to estimate: 



As discussed in section 3T| , the magnitude of the signal depends on the ratio between 
the number of fermions in the fluid and the number of participating nucleons. The 
number of participating nucleons is 0(100), while for having a fluid-like collective 
motion one needs 0(1000) particles, we therefore estimate the required ratio as O(10). 
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We take (j, B , Hi , US ~ 10 MeV~ 0.05/m _1 [| 
The anomaly coefficient is the sum of two terms 



The second term can be neglected because in high energy HIC the chemical energy 
(up) is negligible with respect to the total energy. Also, the axial density is zero at 
to and therefore this term will not effect the initial conditions. 

The quantity that we have to compute is the axial density at freeze-out. This could be 
obtained using the linear approximation described above only if the time of freeze-out 
were in the regime of validity of this linearization. Since the linear approximation 
is not valid for the entire process, one must solve the full set of equations for the 
velocities, energy densities and magnetic fields, and then use them as input for the 
axial current EOM. This task will be left for future analysis. In the following we will 
assume that the general trend (the axial charge being concentrated along the rotation 
axis) will remain the same during the evolution of the system. This assumption may 
be supported by the fact that all the terms in eqs. ( [3,16 ) and ( |3.20|) may be reduced 



due to dissipative effects, but will not change signs. We will approximate the freeze- 
out time by At 10~ 22 sec = 30 fm/c. 

• The boost factor in HIC collisions is taken to be 7 = 100. 

• When integrating over the volume we assume a thickness of order — . 

Before showing the numerical results, let us compare the three QED source terms for 
the axial current: 

'e 2 Z 

CE^Bf ~ C-fFF ~ C7 "~ 

-B 



2 Z\ 1 



Rl J Ro 



~ CfjL 2 d x udtu ~ <V— 5 (3.22) 
QRq 

(where £ is a small factor that takes into account that the relevant regions are close to 
the rotation axis, and not at a distance Ro). We see that the first is larger than the 

second by a factor of ~ 25, and by a factor of 7C {jj^j ~ C " 10 5 than the third. The 
dominant contribution will therefore come from the E-B term. It is worth mentioning that 
although the dominant terms are linearly proportional to the boost factor 7, the volume 
of integration is inversely proportional to 7, and therefore the total axial charge should be 
independent of the collision energy. 

The QCD contribution of the chromo-magnetic field is much more difficult to estimate. 
We will therefore assume that in the deconfined phase, the dominant contribution to the 
colored interaction is the 1-gluon exchange, and that the interaction is similar to the 
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electromagnetic Coulomb interaction, up to a change of coupling constants and group 
theory factors. In this case, the contribution of the chromo-magnetic field to the source 
term of the axial charge is similar in its spatial distribution to the electromagnetic one. 
Assuming that this method of finding the chromo-magnetic field is valid, the ratio between 
the external chromo-magnetic and magnetic field contributions is ^ ~ 100. Finally, note 
that we are not considering topological effects that can induce a change in the total axial 
charge. In the whole discussion, the total axial charge is zero, and we only study the 
implications of its distribution. 

3.4 Numerical Results 

In the following plots we show the numerical results of this analysis for the QED effects. 
As discussed above, we do not have precise values for the various anomaly coefficients and 
external fields and therefore we cannot accurately add the various contributions. However, 
the general features of the effect (axial charge distribution, centrality dependence) are 
similar for all types of contributions, thus combining these results will affect only the 
over-all magnitude. 

In plot | we see the resulting axial charge density as a function of location in the 
plane transverse to the beam direction. We see that the areas of largest charge density are 
located along the axis of angular momentum, and that along the x-axis the charge is zero. 
Plots ||, H and || demonstrate the dependence of the density on the centrality. As expected, 



N participants \p A \ (d^ Term) \p A \ (d^ Term) \p A \ (E^B* Term) 




-4 -2 2 4 -4 -2 2 4 -4 -2 2 4 -4 -2 2 4 

x [fm] x [fm] x [fill] x [fm] 



Figure 3: The number of participants and axial charge density (dark shade means larger absolute 
values) at t — to + At for a midcentral collision (b = Rq). 

the effect is small for central collisions, because of the low angular momentum. 

An important feature of the axial charge distribution is the fact that it is concentrated 
along the rotation axis. In plot ^ we show the second moment of the angular distribution 



/ f dxdvpA(x,v) arctan 2 - 

<^-W j ' (3 ' 23) 
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and its dependence on the centrality (0 is denned as the angle with respect to the x-axis). 
As will be discussed below, this parameter is relevant for the detectability of the proposed 
signal. 
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Figure 4: The axial density for t = to + At, x — for various impact parameters. 
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Figure 5: The axial charge in the upper cone at t = to + At as a function of centrality. The cone is 
centered along the rotation axis, with angular radius of A9 cone = 60°. The centrality is measured 
by the number of participants estimated using the Glauber model. 
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Figure 6: The axial charge in the cone at t = to + At as a function of centrality, divided by 
the number of participants. The cone is centered along the rotation axis, with angular radius of 
A8 cone = 60°. 



3.5 Signal Detection 



The signal discussed in section 3A is proportional to the axial charge in a cone covered by 
the angle Q. According to the results of section ^J, the angles with largest axial charge are 
centered along the angular momentum axis. According to figure ^, most of the axial charge 
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Figure 7: The left plot show the ratio between the axial charge in the belt and in the cone with 
A9 cone = 60° and A9b e i t = 30° (see fig. [l]). The plot on the right shows the second moment of the 
angular distribution defined in [3.23] . 

is located within a cone of angular radius A9 cone ~ 30 — 40° . Plot || also demonstrates that 
matter scattered in the x direction passes through a region of zero axial charge, and should 
have regular production ratios. The scattering into the belt around the equator, with width 
A9b e it can therefore be used to measure ./Va(^(q 4 =o))- The two relevant regions (the one 
with maximal axial charge and the one with zero axial charge) were shown in figure [l|. 

The numerical computation described here can only provide us with the qualitative 
features of the effect under the linear evolution approximation. It does not give us the 
precise axial charge distribution at the moment of freeze-out. This means that it is not 
clear at this stage of the analysis what are the values of A8 cone and AO belt which will give 
the most significant signal. 

Another difficulty in the search for this signature is the fact that unflavored high-spin 
hadrons are very short lived, and decay before they reach the detector. As a solution we 
suggest to focus on narrower high-spin hadrons such as Q~ (see e.g. [|l0|). The angular 
distribution of the production of Q~ can be affected by charge and strangeness effects. In 
order to isolate the spin-dependent effects, it is possible to use E~ as the representative of 
the spin | baryons. This choice of hadrons will decrease the statistics, but will also reduce 
the theoretical uncertainties. J/ip is yet another narrow high-spin hadron which could be 
used as a probe for this signal, but in the case of the charm quark, it is not obvious whether 
the zero mass limit is valid. 

The quantity we therefore propose to measure and compare to the plot is 



N n - (cone) /N n - (belt) _ 
N 3 - (cone) /N s - (belt) 



(3.24) 



as a function of centrality, for various collision energies. As discussed above, this quantity 
is expected to be proportional to ®^ rt , and we therefore expect to find the functional 
behavior seen in plot ||, with no dependence on collision energy. 



4. Discussion 

We presented a proposal for detecting effects of triangle anomalies in heavy ion collisions. 
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It is curious that such subtle quantum effects might be revealed in a collective fluid-like mo- 
tion. Several parts of the presented analysis require careful study. The dynamics of phase 
transitions and the nature of the colored currents in the fluid require a solid theoretical 
understanding. Of particular importance is a precise estimate of the magnitude of the ef- 
fect, which requires a numerical solution to the equations beyond the linear approximation. 
Several important assumptions need clarification: 

• We assumed that in the deconfined phase the gluon field and its contribution to the 
axial current can be estimated assuming a Coulomb interaction. 

• The estimate of the chemical potentials and the assumption that they are coordinate 
independent require a better study. 

• The distribution of the strangeness current can also be centrality dependent. We 
assumed that comparing Q~~ (sss) production with E~ (ssd) production can take this 
effect into account, and isolate the spin dependence, but this needs further study. 

• In the current theoretical framework the freeze-out is successfully modeled by a sharp 
transition from a hydrodynamics description of deconfined QCD fluid with chiral 
symmetry restoration into a kinetic theory of hadron gas. However, if chiral symmetry 
breaking occurs before freeze-out then the enhanced spin alignment discussed here 
might be washed out before the quarks bind into hadrons. The study of this possibility 
requires a better understanding of the phase transition process. 

• The assumption that the fermion number density is proportional to the number of 
participating nucleons requires a stronger theoretical justification. The proportion- 
ality factor is a crucial ingredient in the analysis. 



• In the analysis of section 3.1 the volume over which the integration is performed was 
not specified in details. Since the process of freeze-out is not understood, it is not 
clear what is the relevant volume. We performed the integral over the entire region 
of non-zero axial charge. Taking the integral over smaller regions should increase the 
effect, because the axial charge is concentrated in regions of low participant density. 

• We neglected the contribution of hadronization of gluons in the QCD fluid. 

The last three items can be studied given a full numerical solution to the hydrodynamics 
equations, and a modification of the freeze-out process computation. In the full analysis, 
the number of particles produced during freeze-out is determined by the off-equilibrium 
distribution functions Q 



f l (x^p»)=g i eMP^/T) 



1+ P/^PuT^ 



2T 2 (e+p) 



(4.1) 



in which gi is the degeneracy of the i'th particle species, and r^ v is the dissipative term 
discussed in (|2.2|). In order to take into account the effect of the axial charge we suggest 
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the following modification: 

f*(x^p») =g*exp(p^/T) 



+ 2T2( e + p) a (e+p) 
_ + 2THe + p) +a Z k 9l (e + p) 



(4.2) 



where a quantity marked by * refers to spin excited hadrons (otherwise it refers to a ground 
state hadron), the sums run over all produced particle species, and a is a parameter that 
can be extracted from the data. This phenomenological modification generates the effect 
discussed above, while keeping the total number of produced particles fixed. The results 
can be translated into a cosine expansion of the scattering cross section 

^ oc (1 + 2v 2i cos(20)) (4.3) 
acp 

and the expected result is a suppression of the v^i parameter for spin-excited hadrons. 
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